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The refinement of the Robert-Bonamy (RB) formalism by considering the line coupling for isotropic 
Raman Q lines of linear molecules developed in our previous study [Q. Ma, C. Boulet, and R. H. 
Tipping, J. Chem. Phys. 139, 034305 (2013)] has been extended to infrared P and R lines. In these 
calculations, the main task is to derive diagonal and off-diagonal matrix elements of the Liouville 
operator iS i — S 2 introduced in the formalism. When one considers the line coupling for isotropic 
Raman Q lines where their initial and final rotational quantum numbers are identical, the derivations 
of off-diagonal elements do not require extra correlation functions of the S operator and their Fourier 
transforms except for those used in deriving diagonal elements. In contrast, the derivations for in- 
frared P and R lines become more difficult because they require a lot of new correlation functions 
and their Fourier transforms. By introducing two dimensional correlation functions labeled by two 
tensor ranks and making variable changes to become even functions, the derivations only require the 
latters’ two dimensional Fourier transforms evaluated at two modulation frequencies characterizing 
the averaged energy gap and the frequency detuning between the two coupled transitions. With the 
coordinate representation, it is easy to accurately derive these two dimensional correlation functions. 
Meanwhile, by using the sampling theory one is able to effectively evaluate their two dimensional 
Fourier transforms. Thus, the obstacles in considering the line coupling for P and R lines have been 
overcome. Numerical calculations have been carried out for the half-widths of both the isotropic Ra- 
man Q lines and the infrared P and R lines of C 2 H 2 broadened by N 2 . In comparison with values 
derived from the RB formalism, new calculated values are significantly reduced and become closer 
to measurements. © 2014 AIP Publishing LLC. [http://dx.doi.org/10. 1063/1.4867417] 


I. INTRODUCTION 

Knowledge of molecular spectroscopy is essential for at- 

mospheric sensing, combustion diagnostic, and other scien- 

tific and engineering applications. Among all the spectro- 

scopic parameters of transitions, to accurately determine the 
pressure broadened half-width and the induced shift is the 

most difficult problem and still remains as a big challenge 
for experimentalists and theorists. Because the ambient at- 

mospheric species, temperatures, and pressures are not al- 

ways amenable to laboratory measurements, or because of 

the large number of transitions possible, one often has to 
rely on theoretical calculations based on different line shape 

theories. 1 “ 5 Among them, the formalism developed by Robert 
and Bonamy in 1979 5 is one of the most widely used meth- 
ods, especially for complicated molecules. The main features 
of this well known Robert-Bonamy (RB) formalism consist of 
applying the linked cluster theorem 6 * * * * to evaluate the Liouville 
scattering S operator and introducing an atom-atom model to 
represent short range interactions. 

Mainly due to its importance in practice, a lot of ef- 
forts have been made for years in order to improve the RB 


formalism. 7-10 Descriptions of these refinements and modi- 
fications have been presented in our previous work, 11 they 
are not repeated here again. Despite all efforts to improve 
the RB formalism such that with somehow adjustments of 
potential parameters, the current RB formalism can yield 
reasonably good agreements with measurements, a series of 
recent papers 1213 have demonstrated that for simpler sys- 
tems where results of the close coupling calculations 4 with 
more sophisticated potentials are available, the RB formal- 
ism significantly overestimates the half-widths. Because there 
is no room to make potential adjustments, the large differ- 
ences clearly mean the RB formalism itself contains severe 
weaknesses. 

Recently, we have found that when Robert and Bonamy 
applied the linked cluster theorem to remove the cutoff ap- 
pearing in Anderson’s theory, 1,2 they have assumed matrix 
elements of an exponential function of the Liouville operator 
iS 1 — S 2 can be replaced by values of an ordinary exponen- 
tial function whose arguments are the matrix elements of this 
operator. 1 1 In other words, they have assumed a behavior of 
an exponential of the operator looks like an ordinary func- 
tion. With this replacement, effects from the non-diagonality 
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nature of this operator are ignored. The validity criterion of 
this approximation is that in comparison with its diagonal ma- 
trix elements, the off-diagonal matrix elements of iS i — S2 
are negligible. Usually, this simplification is referred to the 
isolated line approximation, the same name of the approxi- 
mation with which one ignores the off-diagonal matrix ele- 
ments of the relaxation operator W in calculating the spectral 
density F(a>). 

In our previous work , 11 we have scrutinized these two 
simplifying practices in developing the RB formalism and 
have pointed out intrinsic differences between them. Here, we 
briefly outline the differences. These two practices deal with 
different subjects: the non-diagonality (within the linespace) 
of the operator exp(iSi — S2) or of the relaxation operator W. 
Since they cause different consequences, we have proposed 
in Ref. 1 1 to name them differently. The line coupling results 
from off-diagonal matrix elements of exp(iSi — S2), whereas 
the line mixing results from off-diagonal elements of the re- 
laxation matrix W. To neglect the line coupling or to neglect 
the line mixing has a completely different criterion, and the 
former is more stringent than the latter. A failure to distin- 
guish their different criteria by applying the less stringent one 
in both the practices could cause large errors. In many cases 
where the line mixing becomes negligible, effects on calcu- 
lated half-widths and shifts from the line coupling still remain 
important and they have been completely ignored. 

In our previous work, we have developed a refined RB 
formalism applicable for isotropic Raman Q lines of lin- 
ear molecules with which one is able to consider effects on 
calculated half-widths from the line coupling . 1 1 In develop- 
ing this formalism, a manipulating technique plays a crucial 
role to simplify the expression for the off-diagonal matrix 
elements of S2, m iddie such that there are no new correlation 
functions and their Fourier transforms required here except 
for those used in deriving the diagonal matrix elements of 
S2, outer, i, S2, outer, t. and S2, m iddie- Because these correlations and 
Fourier transforms are labeled by two tensor ranks , 10 their 
numbers are only a few. It is worth mentioning that the ma- 
nipulation relies on a condition of &>,y = coff which becomes 
possible only for those coupled lines such that the energy dif- 
ferences between their initial states and between their final 
states are identical (i.e., E? — E, = Ef — Ef). It is obvious 
that for isotropic Raman Q lines of linear molecules, this con- 
dition is exact within the rigid rotor limit and remains well 
applicable when the vibration rotation coupling is taken into 
account. Then, by using this simple expression, we have de- 
rived the off-diagonal elements of S2,middie and subsequently, 
we have successfully considered effects on calculated half- 
widths of isotropic Raman Q lines for N2 broadened by N2 
from the line coupling. The results have clearly demonstrated 
the importance of the line coupling . 11 

Despite this success, there is a weakness remaining in 
this previous study because the simpler expression is only ap- 
plicable for isotropic Raman Q lines. For infrared P and R 
lines where oty,- ^ a>ff , one has to use a general expression 
in which there are a lot of new correlation functions and their 
complex Fourier transforms introduced. Because these new 
functions are not only labeled by the two tensor ranks, but 
also by other parameters, their numbers increase significantly. 


Furthermore, in modern half-width calculations where more 
accurate potential models and the “exact” collisional trajec- 
tory model 7 ' s are adopted, all of these functions have to be 
numerically evaluated. As a result, one faces a challenge in 
deriving off-diagonal matrix elements of S2, middle for P and R 
lines from the general expression. The problem is such serious 
that unless one finds a way to simplify this general expression, 
it becomes very difficult to consider the line coupling for in- 
frared P and R lines at all. Fortunately, we have found that by 
fully exploiting the coordinate representation, one can intro- 
duce two dimensional ( 2 D) correlation functions labeled by 
the same two tensor ranks. With properly changing variables, 
one can make them as even functions for both their new vari- 
ables. Besides, the number of these 2 D symmetric functions 
is very limited. Finally, in terms of the latters’ 2 D Fourier 
transforms, one is able to obtain a simpler expression for the 
off-diagonal elements of S2, m iddie- 

Besides obtaining the simpler expression, one enjoys 
other benefits by easily discovering some hiding features of 
the off-diagonal elements. For example, because the 2 D cor- 
relation functions are even functions, their 2 D Fourier trans- 
forms become real. Then, it becomes obvious that the off- 
diagonal matrix elements of S2, middle are real too. In addition, 
the averaged energy gap and the frequency detuning between 
the two coupled lines appear as the two modulation frequen- 
cies of these 2 D Fourier transforms in the new expression. 
Based on this feature, one can conclude that it is these two 
quantities that characterize the coupling strength between the 
two coupled lines. Armed with this knowledge, to predict 
whether the coupling between the two coupled lines is im- 
portant or not becomes easier and more accurate. 

Of course, the success in applying the simpler expres- 
sion relies on how to effectively and accurately evaluate these 
2 D Fourier transforms in practical calculations. In our pre- 
vious works, based on the sampling theory we have devel- 
oped effective tools to derive one dimensional Fourier trans- 
forms from one dimensional correlation functions. By com- 
bining the techniques with other approaches available in the 
signal theory , 14 we are able to successfully derive these 2 D 
Fourier transforms. After the latter are available, there are no 
more obstacles left in considering the line coupling for P and 
R lines. Numerical calculations of the half-widths for the in- 
frared P and R lines of C2H2 broadened by N2 based a new 
potential model 1 '’ together with those for the isotropic Raman 
Q lines have been carried out. The calculated results clearly 
demonstrate that effects from the line coupling are important 
in both these cases. In comparison with values derived from 
the RB formalism, new calculated values are significantly re- 
duced and become closer to measurements . 16 

In Sec. II, we show how to introduce the 2 D correlation 
functions and to change them to become even functions. By 
exploiting the latters’ symmetry properties, a simpler expres- 
sion for the off-diagonal matrix elements of S2, middle can be 
obtained in terms of their 2 D Fourier transforms. With this ex- 
pression, to calculate the off-diagonal matrix elements for P 
and R lines becomes more tractable in practical calculations. 
In Sec. Ill, we apply the new formula to the C2H2-N2 sys- 
tem and present calculated half-widths of the infrared R lines. 
In Sec. IV, we present discussions and conclusions. With 
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detailed formula derivations and numerical sample demon- 
strations presented in Secs. II and III, we believe our refine- 
ment of the RB formalism makes another progress in the right 
direction. 


II. THEORY 

A. Matrices of S 2 , outer, i; outer, fj 3nd S2, middle 

in the coordinate representation 

It has been known for years, when one calculates ma- 
trices of S 2 , outer, i, s 2 , outer, f, and Saddle based on potentials 
containing site-site model components, one could encounter 
convergence problems. 17 As demonstrated in our previous 


work, 10 a solution is to use the coordinate representation by 
choosing the orientations of the pair of molecules as the basis 
set in Hilbert space; i.e., |<5.(f2 fl — £2 no ,)} ® |.<5(f2/ a — i, a )), 

where £l m and Q| w represent orientations of the absorber and 
bath molecules specified by a, respectively. In contrast, with 
the standard representation, the basis set is constructed from 
<g> ].i 2 m 2 ), the product of the states of two interacting 
molecules. The advantage in choosing the coordinate repre- 
sentation results from a fact that with this representation, the 
potential becomes a diagonal operator and its matrix elements 
become multi-dimensional integrations whose angular parts 
can be evaluated analytically. 

With the coordinate representation, one is able to write 
the matrix elements of S 2 i 0 u tei-,i as 11 


OO t 

Sxou/erM) = 1 2 f dt f dt' « if' JMj\{L\(t)L\(t'))\if, JMj 

/ OO p t 

dt / dt ' E Pi2 

-OO j OO • m :/ 


», 


outer, i 


-OO — OO 




i2tn2 nti i"m " i 2 m 2 


n 2 (2ji + 1) 

x J d£ 2 a J d£2fi{i f mji2m2\(x)V a (R(t))(a\i"m / j 'i' 0 m, 

x {i"m"i' 2 m' 2 \f})Vp(R{t')){f}\imii 2 ni 2 ), 




(1) 

where can" = [£^(0 — Ef a Xi")\Hi, a)i 2 q = [£ (fc) (i 2 ) — and la) is a short notation for |<5(f2 fl — <S> 

— Qba)) and V a (R(t)) represents the potential evaluated at this specified orientation labeled by a. Similarly, one is able to obtain 
the expression for S\ { iu l / er f(r c ) which is the same as Eq. (1 ) except for making exchanges between the quantum numbers i and 


f' , between i' and f, and making switches between t and t'. Meanwhile, one can obtain an expression for S' 


i'f'Jf 

2 , middle 


( r c ) as 1 


Q 1 ' f'df ( ) 

^Iwiddlevc) ~ *0 


00 1 

= p / * / <"’ K< (f ' 


\(t'))\if, JMj middle 


— OO — OO 


OO OO 


m 1 


1 m 1' m' 1 m \ 




— OO — OO 


i 2 m 2 i' 2 m ’ 2 (m) 

x (-1 m\ - m' f Mj)C(if J, m, - m f Mj ) 

x J d£l a J d£ip(i'm' i i2in2\oi)V c ,(R(t))(ci\imji2m2) 

x {fm f i 2 m 2 \P)Vfi(R(t')){P\f'in' f i' 2 m' 2 ). 


(2) 


It is worth mentioning that as one considers the line coupling 
occurring among lines within the same bands, expressions for 
the matrix elements of S 2 : 0 ute r,i and S 2 oute r,f can be simpli- 
fied further. It is obvious that among lines within the same 
bands, both their initial and final vibrational quantum num- 
bers are identical. This implies that one can replace Sjjj. by 

Si’i in the expression for S‘ 2 { m l / er ,■ (f c ) and replace 8 Jf f by 8g 


in ^2 mitlr /( r c)- As a result, the matrices of both these two 
terms are diagonal. But, this simplicity does not apply for the 
matrix of S 2 . middle which remains as an off-diagonal matrix 
even within the same bands. 

For systems consisting of two linear molecules, formulas 
used to evaluate the diagonal matrix elements of the S 2 iOUter ,i, 
S 2 , outer, f, and S 2jln iddie terms are well known. 10 With respect to 
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expressions for the off-diagonal matrix elements of Sa.middie, a 
simpler formula applicable for isotropic Raman Q lines with 
(Oft = (Off and a more general one applicable for lines in the 
P and R branches with con ^ a>ff are also available. 1 1 How- 
ever, the more general formula is somehow awkward in practi- 
cal applications. It becomes desirable to find a more tractable 
expression. 

B. Expression for matrix elements of S2, middle 
in the P and R lines 

Usually, potential models are given in terms of a spherical 
tensor expansion as 10 


V(R(t))= U(L x L 2 L-,R{t)) 

L\L 2 L 

X ^2 C{L\L 2 L, mim 2 m)Y Limi (Q. a ) 

m\m 2 m 

x Y L2m2 (Q h )Yl m (w(t)). (3) 

By inserting this expression for the potential into Eq. (2) and 
analytically carrying out integrations over the orientations of 
the two molecules and completing summations over the mag- 
netic quantum numbers as more as possible, one is able to 
obtain a new expression for the matrix elements of S2, middle 
as 


SZ£JLe<rc) = + l)(2j' + 1)(2 ji + 1)(2 j f + 1) 

X £(-l) 1+/+Ll W(j Wfjijf’ JLACUuiLu Q00)CUfjfLi, 000) 


L\L 2 


X J2 (2l i + !)(2*2 + l)pi 2 C 2 {i 2 i' 2 L 2 , 000) 


1 


l 2 l 2 


16jt 2 ^ 2 (2L 1 + 1)(2L 2 + 1) 


OO OO 

J2(-l) L>+L > +L J2 f dt f dt’e i{a>i ' i+m ‘ W , e - i(m <' s ^ y 

T Ayf J J 


M 


— OO — OO 


x U {L\L 2 L\ R(t))U (L\L 2 L\ R(t'))Y* M {6 a {t), (j> a (.t))Y LM (e p {t'), <^(f')). 
Then, one can introduce 2D correlation functions defined by 

1 


G LlLl (t, t') = 


l6n 2 ti 2 (2Li + 1)(2L 2 + 1) 


D-d 


L\+L 2 +L 


x U(L x L 2 L ■ R(t))U (L\L 2 L\ R(t'))J2Y* LM (^(t), 4> a {t))Y LM (0p{t\ </>p(t')) 

M 

1 

(4^) 3 ^ 2 (2Li + 1)(2L 2 + 1) 


^(_1)L i +l 2 +L(2 L+ 


x U (L\L 2 L\ R(t))U {L\L 2 L\ R(t'))P L ( cos©,.,-). 


where 


(4) 


(5) 


cos©, = cos 0 a (t)cos6p(t') + sin# a (/)sin 9p(t'). (6) 

In deriving the above expression, we have assumed that collisional trajectories lie in a half part of the x and z plane with positive 
x values. As a result, one can simply set (p( t) = (p( t') = 0. Then, in terms of these 2D correlation functions, the off-diagonal 
matrix elements of Si. middle can be expressed as 

sZfL,e(rc) = (-iy f+j ^(2j; + 1)(2 j'f + 1X2./, + 1)(2 jf + 1) 

X £(-l) l+1+Llw U'ii'fiiif, JLOCUijlLu 000)CUfjfLu 000) 

L\L 2 

x J2 (2l 2 + 1)(2 i 2 + \) Pil C 2 {i 2 i' 2 L 2 , 000) 

* 2*2 


OO OO 

J I t'). 


(7) 


At this stage, unless con — (Off, our previous manipulation technique 10 - 11 used in deriving the diagonal matrix elements 
of S 2 , 0 uter,i, S 2>oute r,f, and S 2>m iddie and also in deriving the off-diagonal matrix elements of S 2, middle for isotropic Raman lines 
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ill the Q branch is not applicable. This implies that for lines in the P and R branches, one has to find other approaches. 
Let us focus attention on the integration part of Eq. (7). By changing its integration variables t and t' to the variables of r 
= t — t' and r' = (l/2)(t + t'), it can be rewritten as 

OO OO 
— OO — OO 


-II 


drdr'e 


i ( o > ,/,■ +01^ )(t '+!)-i (fOfif+eo^ , 2 )(t'~ \ ) 


! 2 ‘ 2 ' 



\) 


-II 


dTdT'e , y-^ +a ^) x+,ia,t, - co f' l ' )r G LlL 2 (r, r'). 


( 8 ) 


where we have introduced the new 2D symmetric correlation 
functions defined by 

Gi,L 2 (r, r') = G LiLi (x' + x' - . (9) 

It is worth mentioning that in terms of new variables r and 
r\ the original ones are given by t = r' + (l/2)r and t' = r' 
— (l/2)r, respectively. 

C. Symmetry properties of the 2D correlation 
functions 

It turns out that by exploiting symmetry properties of 
these new 2D correlation functions, one is able to sim- 
plify Eq. (8) further. With respect to the original correlation 
functions G^lAT t'), they have symmetries of GnL 2 (t, t') 
= G Ll L 2 (t', t) and G LlLl (t, t') = G LlL2 (-t, -t'). The latter 
represents the fact that collision processes are time-reversal 
invariant. But, t') are neither even functions of t nor 

even functions of t'. Unfortunately, a failure in pertaining to 
the evenness for both their two variables t and t' would block 
their usefulness. 

On the other hand, one can conclude that the new cor- 
relation functions G^/^r, r') are invariable for an operation 
of r -o- — r, an operation of r' -o- — r', and an operation for 
both r -o- — r and r' -o- — r'. The first two properties mean 
they become even functions both over r and r'. The last one 
represents the fact that collisional processes are time-reversal 
invariant. The above conclusions are based on the symmetry 
properties of the functions of G/, 1 L 2 (r' + §, t' — |), 

Gl ^ ( r ' + ^ ,r ' - D 

= G LlLl (- r ' - X ~) 

= G LlL2 (-r'- T -,-r' + l). (10) 

Meanwhile, in contrast with G [ t L 2 (t, t'), G^^Cr, r') do not 
have a symmetry between r and r'. In other words, they vary 


differently as r and r' vary. It is worth mentioning that all 
these 2D correlation functions are associated with specified 
collisional trajectories. Usually, the latter are labeled by the 
closest distance r c (or the impact parameter b). As a result, 
they are also functions of r c . For simplicity, this dependence 
may not explicitly be presented. 

As example, we consider the C 2 H 2 -N 2 pair in the present 
study. If one assumes no bending modes of the C 2 H 2 molecule 
are excited and only considers its stretching bands, this 
molecule remains linear. Thus, all the formulas developed 
above are applicable for this system. Besides, a new potential 
model containing 85 spherical components labeled by com- 
binations of even values of the tensor ranks Li, L 2 , and L is 
available recently in the literature. 1 ^ Based on this potential 
model, one can calculate the original 2D correlation functions 
G /, 2 U , t') associated with a specified selection of Lj and L 2 
and the newly introduced ones Gl 1 l 2 (t, x '). In practical cal- 
culations, instead of t and t' people prefer to use dimension- 
less variables z = vt/r c and z! = vf/r c , where v is the averaged 
velocity for the “exact” trajectory model and is the “apparent” 
velocity for the “parabolic” trajectory model. 5 We follow this 
common custom here. 

Thus, there are two sets of variables used in the present 
study. One set consists of variables with dimensions that ap- 
pear in formulas and another set consists of dimensionless 
variables which appear in plots. Within the first set, there are 
original t and t' and new introduced r and r'. Within the di- 
mensionless set, there are z and z' corresponding to t and 
t', and there are u = vr/r, and v = vx' / r c corresponding to 
r and r'. It is obvious that u = z — z' and v = (z + z')/2. 
For simplifying notations, we use the same function sym- 
bol for G LlLl (t, t') and G^ jl 2 (z, z'), and the same symbol for 
Gl!L 2 (t, t') and G/,jL 2 (m, v ). In order to avoid confusion, we 
list notations used in the present study in Table I. 

In order to show profiles of Gl 1 l 2 (z, z') and Gl 1 l 2 (m, v ), 
we consider their most important components associated with 
L! = 2 and L 2 = 0 and present their three dimensional profiles 
obtained at r c = 4.5 A in Figs. 1 and 2, respectively. To choose 
r c = 4.5 A as an example mainly results from a fact that the 
profiles at this region contain more structures and their profile 
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TABLE I. Notations used in the present study. 


Variables with dimension 


Dimensionless variables 


Correlation functions 

Fourier transforms 

Correlation functions 

Fourier transforms 

/') 

GijLj (r, r'), where r = t — t' and r' = (t + t')/2 

®LiL 2 (<», <*>') 

G Li l 2 (z,z') 

G Li l 2 (u, v), where u = z — z' and v = (z + z')/2 

F Ll L 2 (k,k') 


differences can be more clearly exhibited. As shown in Fig. 1, 
Gioiz, z') is neither even over z nor even over z! . In contrast, 
G 20 (u, v ) becomes an even function over both u and v. 

D. The 2D Fourier transforms of the 2D correlation 
functions 

Because r, r') become even functions of r and r', 

their Fourier transforms are real. One can introduce their 2D 


Fourier transforms defined by 

OO OO 

F LlL2 (<w, a/) = ^ J J dr dr' e lmz e ,a> ' z ' G LlLl (j , r'). 

— OO — OO 

( 11 ) 

Then, in terms of these Fourier transforms, one is able to 
rewrite Eq. (8) as 


OO OO 


J /» OO 

J — 00 J — c 


dzdz e 


— +".u 2 J r +‘ (®/t > r 


-OO «/ — OO 
OO OO 


/ OO /»C 
-OO t/ — c 


dcodcoe la>z e ,a>z F LlL2 (co, co ) 


OO OO 

^ J J dcodw'W LlL2 (a),(o') J J 


dzdz e 


( (w.'/j+Wf//) \ 

2 7 7 +<W y / ,- 2 - T + f (<W /f -<W yV,-/ ~(o')Z 


— OO — OO 
OO OO 


~ 27T / / — - W/J 12 ) - <W/i + &>/'/') 


= 27rF Lli . 2 


(Oi’i + C Off 


+ «d; 2 , w fi' — <W/i I • 


( 12 ) 


Finally, in terms of these 2D Fourier transforms F l,l 2 (®, co'), the off-diagonal matrix elements of S 2 , m iddie can be simplified as 

<£&;,frc) = 27 r (— \y t+ j'f y (2j- + d(2 r f + mj, + mjf + d 

x J2 ( - V) ' +J+L ' W( J‘jfJ‘jf’ JL^CUijlLuOmCijfjfLuOOO) 


L1L2 


x £(2i 2 + l)(2i _ 2 + l)pi 2 C 2 ( 121 ^ 2 , 000)F£ji 2 


I 2 ! 2 


"i'i + «/'/ 


“b tU(' l2 , (Of'i ' (T)fi J ■ 


(13) 


In comparison with the general expression derived in our 
previous study 11 and other formulas available in literature, 20 
Eq. (13) has several advantages. First of all, this expression 
is more concise. Besides, because the 2D Fourier transforms 
Fl 1 l 2 ( co ’ c ° ) are rea l’ one can easily draw a conclusion that 
the off-diagonal elements of S2, middle are real too. In con- 
trast, with other formulas, to verify the above conclusion 
is not such straightforward because in general, they contain 
complex functions. However, the most advantage in usage of 
Eq. (13) is that the number of functions F LlLl (co, co') is very 
limited. As example, for the system of C 2 H 2 -N 2 , by using 


Eq. (13) one only needs to evaluate six 2D Fourier trans- 
forms labeled by Li and L 2 as (20), (22), (24), (40), (42), 
and (44). In practical calculations, in order to simplify fur- 
ther, one can omit three minor ones with (24), (42), and (44). 
It is well known that as long as the “exact” trajectory model is 
adopted in calculations, all the functions introduced such as 
those resonance functions appearing in literature have to be 
numerically evaluated. If one uses other formulas to calculate 
the off-diagonal matrix elements of S 2 , middle* the number of 
these functions would increase dramatically. In addition, they 
are complex functions. 
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FIG. 1 . The 2D correlation function of Giofa <■ ) (in ps" 2 ) at T = 296 K and 
r c = 4.5 A for a molecular pair of C 2 H 2 -N 2 . The calculation is based on the 
new potential model 13 and the “exact” trajectory model. 

It is interesting to note that there is a simple relation be- 
tween the 2D correlation functions and their one dimensional 
(ID) partners introduced previously 10 in calculating the diag- 
onal matrix elements of S 2 , 0 uter,i, S 2 , outer, f, and S 2 , m iddie- With 
the current notations, the ID correlation functions introduced 
previously are defined as 

OO 

G Li lA t)= J dz'G LiLl {T,x'). (14) 

— OO 

Similarly, there is also a simple relation between the 2D 
Fourier transforms F LlLl (a>, co') and the ID Fourier trans- 
forms Fi 1 l 2 (®) which is defined by 

OO 

F Ll L 2 M = -L j dxe iar G LlLl {x). (15) 

— OO 

Because co corresponds to r and the latter represents the gap 
between t and t', it becomes clear that it is how the original 
2D correlation functions G / /, 2 (f , t') vary with t — t' that mat- 
ters at the last stage of calculations for the diagonal matrix 
elements of S 2 ,outer,i> S 2 , outer, f? and S 2 , middle- 
Finally, it is easy to find out that 

OO OO 

F Li l 2 (co, 0) = ^ J dxe l(OX J dz'G LlLl (z, x') 

— OO — OO 

OO 

= ^ J dxe ia>r G LlL2 (r) 

— OO 

= /— F L{Ll (co) . (16) 

V Z.71 

As a result, it becomes obvious that the expression for the off- 
diagonal matrix elements of S 2 , m iddie is also applicable for its 
diagonal matrix elements. 


FIG. 2. The same as Fig. 1 except for the 2D correlation function of 
G 20 

E. Tools used in numerical calculations 

In practice, the success in applying the new method to 
calculate these off-diagonal matrix elements of S 2 , middle relies 
on how to effectively and accurately evaluate these 2D Fourier 
transforms. Thanks to our experiences in calculating pressure 
broadened half-widths and induced shifts, we have accumu- 
lated enough tools to overcome this challenge. First of all, 
by introducing the coordinate representation, the main calcu- 
lation task is reduced to calculate the 2D correlation func- 
tions and their 2D Fourier transforms. As mentioned above, 
the number of these functions is very limited. In addition, as 
shown in Eq. (5), no matter how complicated the potential 
model used would be, there is no any difficulty to evaluate 
these 2D correlation functions accurately. 

Second, in our previous studies we have developed an 
effective approach to derive ID Fourier transforms. This 
method is based on the sampling theory. 21 The basic idea is 
by sampling, one converts a function of interest to a sequence 
in the time domain, calculates its discrete Fourier transform 
with the fast Fourier transforms (FFT), 22 and then relates this 
transformed sequence in the frequency domain to the Fourier 
transform of the original function. 

When we try to extend our study on the line coupling 
from Raman Q lines into infrared P and R lines, we have faced 
a challenge how to evaluate a 2D Fourier transform from a 2D 
correlation function. Fortunately, in other scientific and engi- 
neering fields such as the signal theory 14 to deal with data pro- 
cessing and analysis, to consider the 2D Fourier transforms of 
the 2D functions is a familiar subject. There are many books 
and papers in literature and there are several tools available 
there with which one can compute the 2D Fourier transforms. 
For example, it has been known that the 2D Fourier transform 
is linearly separable: the Fourier transform of a 2D image 
is the Fourier transform of the rows followed by the Fourier 
transforms of the resulting columns (or vice versa). Based on 
this fact, the 2D Fourier transform can be carried out by ID 
Fourier transforming all the rows of the 2D signal and then ID 
Fourier transforming all the columns of the resulting matrix. 
In the present study, because we have very effective programs 
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FIG. 3. The 2D Fourier transform F 2 o(£, k') (in ps -2 ) at T = 296 K for a 
molecular pair of C 2 H 2 -N 2 . The calculation is based on the “exact” trajectory 
model with r c = 4.0 A. 

to carry out the ID Fourier transform, we have adopted this 
straightforward way to compute the 2D Fourier transform and 
have obtained satisfactory results. 

Based on the new potential model, 15 one can find the 
minimum of the closest distance for the “exact” trajectories 
at T = 296 K is r cm ; n = 3.7546 A. Because nearly head- 
on collisions play a dominant role in the line coupling, we 
select a trajectory with r c = 4.0 A as example. Among all 
the components, we consider the most important one with Li 
= 2 and L 2 = 0. Again, we prefer to plot functions with di- 
mensionless arguments used in our numerical calculations. In 
this case, the function is FioiA, k f ) which is the Fourier trans- 
form of G 20 (u, v) . Thus, we present the three dimensional 


plot of F 2 o(A, k') obtained at r c = 4.0 A in Fig. 3. This plot 
provides a whole picture about how its magnitudes vary with 
its two variables k and k' at this specified trajectory. 

In order to more clearly provide some numerical mea- 
sures on magnitudes of IF 20 ( A , k'), we present another plot to 
show its profiles as functions of k or k' with three fixed k' or 
k values of 0.0, 8.0, and 16.0 in Fig. 4. By combining Figs. 3 
and 4, one can find that in terms of the coordinates of k and 
k', there are four peaks with their heights at (0, 10), (0, —10), 
(5, 0), and (—5, 0). The first two peaks are symmetrically 
located at k' axis and the last two are symmetrically at k 
axis. In addition, the former’s height is almost identical to the 
latter’s height. But, shapes of the peaks in different pairs are 
quite different. In general, the peaks in the first pair distribute 
more widely than those in the second pair. As shown in 
Fig. 4, roughly speaking after k > 15, magnitudes of 
F 2 o(k, k') would decrease more than one order from its 
maximum. On the other hand, its magnitudes could still 
remain significant or non-negligible as long as k' < 30. We 
will be back to this subject later. 

III. CALCULATIONS OF HALF-WIDTHS FOR C 2 H 2 
LINES BROADENED BY N 2 

A. Half-widths of isotropic Raman Q lines 

First of all, we present our calculated half-widths for 
isotropic Raman Q lines of C 2 H 2 broadened by N 2 obtained 
from including the line coupling. The calculation is similar to 
that reported in our previous work for N 2 lines broadened by 
N 2. 11 Here, we briefly outline some of main factors to deter- 
mine effects from the line coupling. First of all, let us look 
at the dashed-dotted black curve with k' = 0 in Fig. 4. Af- 
ter reaching its maximum at k = 5, magnitudes of this curve 



FIG. 4. F 2 o(A k') (in ps 2 ) at r c = 4.0 A and T = 296 K with three fixed k' values as functions of k, plotted by three dashed-dotted curves. Meanwhile, 
F 2 o(A k') as function of k' with three fixed k values are plotted by three solid curves. The black, red, and green colors correspond to the three fixed values of 
0.0, 8.0, and 16.0. 
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FIG. 5. Comparison of calculated half-widths at 296 K for the isotropic Raman Q lines of C 2 H 2 broadened by N 2 . Theoretically calculated results from the 
modified RB formalism without and with the line coupling are plotted by + (red) and (green), respectively. Meanwhile, measured values at 298 K 23 are 
given by x . 
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decrease as k increases, but still remain significant within a 
range of k < 14. This curve represents the profile of F 20 ( k , 0) 
and the latter is associated with F->o ( a> , 0) (= -J= F 20 (a>)). For 

the specified conditions of the plot (i.e., r c = 4.0 A and T 
= 296 K), the value convert from co (in cm -1 ) to k (dimen- 
sionless) is k ?=» 0. 1 1 x co. Thus, we know that as long as 
the first argument of F 20 0) in Eq. (13) is less 

than 127 cm' 1 , one has to consider the line coupling. On the 
other hand, the rotational constant B of the C 2 H 2 molecule is 
only around 1.2 cm -1 . This implies that for the most nearby 
coupled lines, their initial (or final) energy gaps together with 


u>i' 2 i 2 are well within this argument range to make contribu- 
tions to the off-diagonal matrix elements of S 2 , middle- 

Based on these facts, one can expect that effects on cal- 
culated half-widths of C 2 H 2 immersed in N 2 bath would be 
significant. In Figs. 5 and 6, we present our calculated values 
by considering the line coupling at T = 296 K and T = 150 K, 
respectively, together with those derived from the modified 
RB formalism without taking into account the line coupling. 
As comparisons, recent measured data 2 ' at T = 298 K and 
T = 150 K are also plotted in these figures. As shown in 
Figs. 5 and 6, by taking into account the line coupling, the new 


1 — ■ — 1 — ' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ' — 1 — ' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ' — 1 — ' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ' — 1 — ' — r 


T = 150 K 



FIG. 6. The same as Fig. 5 except that calculations and measurements 23 are earned out at 150 K. 
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TABLE II. Energy and frequency values for coupled P and R lines. 


Lines 

P(j— 2) 

j - 3 <-j - 2 

R(j-2) 

j — 1 ■*- j — 2 

P0) 

j - 1 «-j 

RO) 

j+i 

P0+2) 

i+1 ^j + 2 

RO + 2) 
j + 3 j +2 

Ef and Ej (in B) 

(j— 3)0— 2), 

O-Bj. 

0-i)j. 

0 + DO + 2), 

0 + DO + 2), 

0 + 3)0 + 4), 


0-2)0- 1) 

0-2)0- 1) 

jO + 1) 

iO + i) 

0 + 2)0 + 3) 

0 + 2)0 + 3) 

(Oft (inB) 

-20-2) 

20-1) 

— 2j 

20 + 1) 

-20 + 2) 

20 + 3) 


calculated values at 296 K and 150 K are reduced by about 
13% and 14%, respectively, and they become closer to the 
measured values. Thus, one can conclude that it is abso- 
lutely necessary to consider the line coupling in calculating 
N 2 broadened half-widths of isotropic Raman Q lines of the 
C 2 H 2 molecule. 


B. Half-widths of infrared P and R lines 

In this section, by considering infrared P and R lines of 
C 2 H 2 immersed in the N 2 bath, we present some things new. 
We would like to mention that we have adjusted our practi- 
cal calculations to match the potential model adopted. 15 First 
of all, we provide some facts. It has been known that for 
lines within the same bands, the S 2 , ou ter,i and S 2 , ou ter,f terms 
are diagonal and their diagonal matrix elements are complex. 
Usually, ImS2,outer,i and ImS2, ou ter,f are one order smaller than 
R e S2,outer,i and ReS2.0uter.f- Meanwhile, the S2, middle term is al- 
ways off-diagonal, but all its matrix elements are real. Finally, 
the iS 1 term is imaginary and it is the main component of ImS. 
If the potential models do not have vibrational dependences, 
iS 1 becomes zero. 

Because the current potential model 15 does not have vi- 
brational dependences, we are not able to properly evaluate 
the iSi term at all. This deficiency would cause a large un- 
certainty in evaluating the whole ImS. Based on this and the 
other fact mentioned above, we have decided to completely 
ignore the whole ImS (=iSi + ImS 2 . 0 uter.i + ImS 2 , outer, f) in 
our practical calculations. 

Before presenting our final results, we would like to pro- 
vide estimation of how important the effects from the line 
coupling would be. The effects are mainly determined by that 
in comparison with the whole diagonal matrix elements of 
ReS 2 , whether the off-diagonal matrix elements of S 2 , m iddie 
given by Eq. (13) are significant or not. By analyzing the pro- 
files of 2D Fourier transforms F^/,, (co, co') and considering 
their two arguments of + Wi ,^ an d &>/y — cof appear- 

ing in Eq. (13), one can find the answer. First of all, because 
there is no Q branch in stretching bands, one only needs to 
consider P and R lines. Second, two coupled lines must share 
the same parity for their initial angular quantum numbers. The 
same conclusion is also true for their final quantum numbers. 


Based on the evenness of the tensor rank Li, the whole line 
space constructed by all the P and R lines can be divided 
into two independent sub-spaces. One consists of R(0), P(2), 
R(2), P(4), R(4), . . . and another consists of P(l), R(l), P(3), 
P(3), and so on. Then, lines only in the same groups can be 
coupled. We list general expressions for their initial and fi- 
nal rotational energies together with their frequency exclud- 
ing vibrational band center for lines in the same groups in 
Table II. With Table II, one can find corresponding averaged 
energy gaps (= av, + a>/ ' / ) anc | frequency gaps (=oof{ — cofi) 
between coupled lines. Besides, the general expressions for 
these two quantities are presented in Table III. 

As shown in Table III, the frequency gaps between two 
coupled nearby R lines are constant (i.e., ±4.8 cm -1 ). In con- 
trast, the frequency gaps between two coupled nearby P and 
R lines are ±2(2j ±1) and ±2(2j ± 3) which increase very 
quickly as j increases. We note that after making the variable 
exchange, this argument appears as k' in the 2D Fourier trans- 
forms k'). For the specified trajectory with r c = 4.0 

A, co — 4.8 cm -1 corresponds to k = 0.53 which is pretty 
small. Then, roughly speaking, contributions to the line cou- 
pling between two nearby R lines are mainly determined by 
profiles of Fl,l 2 (k, 0) . Among them, the profile of their ma- 
jor component F 20 (k, 0) is given by the dashed-dotted black 
line in Fig. 4. In contrast, contributions to the line coupling 
between two nearby P and R lines are from other F LlLl (k, k') 
with k' > 0 whose magnitudes are smaller than F LlLl ( k , 0) . 
The larger the j value is, the larger the k' is and the smaller 
the contributions are. Therefore, one can expect that in each 
of these two sub-spaces, especially for lines with large j val- 
ues, the line couplings mainly occur among lines within the 
same branches. In other words, in comparison with the cou- 
plings within the same branches, except for lines with small j 
values, the line couplings between the P and R branches are 
not important. Our numerical calculations shown later have 
verified this expectation. This result was also previously ob- 
tained from fully quantum calculations of line mixing effect 
in CO-He. 18 

For each of selected 600 values of r c ranging from r c mln to 
40.0 A, we have calculated all the diagonal and off-diagonal 
matrix elements of S 2 (r c ) in the sub-space constructed by cou- 
pled lines of R(0), P(2), R( 2 ), . . . . , P(40), R(40), and P(42). 
The size of these matrices is 42 x 42 which is large enough. 


TABLE III. Averaged energy gaps and frequency gaps between coupled nearby lines. 


Coupled lines 

R0 _ 2) ~ RO) 

P(j) ~ RO) 

P0 + 2) ~ RO) 

RO + 2) ~ R0) 

Energy gap (in B) 

±4j 

±(2j + 1) 

±(2j + 3) 

±40 + 2) 

Frequency gap (in B) 

±4 

±2(2j + 1) 

±2(2j + 3) 

±4 
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but remains tractable. Because we have set 4 as the upper limit 
of Li, these S 2 (r c ) become five diagonal matrices. In addi- 
tion, they are real and asymmetric. With the same method de- 
veloped in our previous work , 11 we can find the matrices of 
exp[— S 2 (r c )]. After the latter are available, by using the ex- 
pression for matrix elements of the relaxation matrix W given 
by 


+oo 



x {Sniff - « i'f\e~ s ^\if »}, (17) 

one can easily derive all the matrix elements of the relax- 
ation operator W whose diagonal matrix elements correspond 
to the calculated half-widths with including effects from the 


line coupling. The results are presented in Matrix 1 where 
due to lack of room to print all its contents, some rows and 
columns are removed. As shown in Matrix 1, the W matrix 
has large and positive diagonal elements representing the cor- 
responding half-widths. With respect to its off-diagonal el- 
ements, those corresponding to the line mixing within the 
same branches (i.e., the P-P and R-R mixings) are negative 
and that associated with the line mixing between two different 
branches (i.e., the P-R and R-P mixings) are positive. In other 
words, in comparison with the P-P and R-R mixing terms, 
the P-R and R-P mixing terms have an opposite sign, as also 
previously observed for CO-He. ,s In addition, the former’s 
magnitudes are significantly larger than the latter’s. This fea- 
ture results from the fact that the line couplings within the 
same branches are more important than that happening be- 
tween different branches. 


Matrix 1 : A 42 x 42 matrix of the relaxation operator W in the P and R line space 


122.064 

7.265 

-10.247 

2.984 

-3.739 

2.072 

-2.454 

1.583 

-1.790 . 

. -0.003 

0.001 

7.216 

115.901 

4.528 

-12.931 

2.876 

-5.121 

2.213 

-3.391 

1.751 . 

. 0.003 

- 0.002 

-10.077 

4.482 

112.361 

3.528 

-13.935 

2.645 

-5.969 

2.128 

-3.987 . 

. -0.005 

0.002 

2.896 

-12.614 

3.483 

109.017 

2.982 

-14.047 

2.404 

-6.467 

1.999 . 

. 0.004 

-0.004 

-3.562 

2.773 

-13.452 

2.941 

105.603 

2.607 

-13.935 

2.192 

-6.762 . 

. -0.008 

0.002 

1.951 

-4.798 

2.541 

-13.362 

2.570 

102.418 

2.323 

-13.926 

2.020 . 

. 0.004 

-0.006 

-2.247 

2.072 

-5.488 

2.306 

-12.998 

2.290 

99.803 

2.103 

-14.191 . 

. - 0.011 

0.003 

1.443 

-3.036 

1.988 

-5.816 

2.102 

-12.689 

2.074 

97.896 

1.933 . 

. 0.004 

-0.008 

-1.577 

1.595 

-3.486 

1.870 

-5.934 

1.939 

-12.606 

1.907 

96.603 . 

. -0.016 

0.003 

- 0.001 

0.001 

- 0.002 

0.001 

-0.003 

0.001 

-0.004 

0.002 

-0.005 . 

. 46.968 

0.037 

0.000 

- 0.000 

0.001 

- 0.001 

0.001 

- 0.002 

0.001 

-0.003 

0.001 . 

. 0.034 

44.758 


Similarly, we have calculated the matrix elements of W in 
the sub-space constructed by other coupled lines of P(l), R(l), 
P(3), . . . , R(39), P(41), and R(41). It turns out that because the 
vibrational dependences have not been taken into account, this 
W matrix is identical to the previous one. As a result, one can 
identify the half-widths not only for all the even R lines, but 
also for all odd R lines from Matrix 1 alone because the half- 
width of 2j < — 2j — 1 (i.e., R(2j — 1)) is the same as 2j — 1 
4 — 2j (i.e., P(2j)) with the rigid rotor approximation . 11 

After all these values are available, we present the cal- 
culated half-widths of the R lines obtained from considering 
all possible line couplings within all P and R lines together 
with those derived from the modified RB formalism without 
the line coupling in Fig. 7. For comparison, some measured 
values 16 are also plotted. As shown in the figure, values ob- 
tained from the RB formalism significantly overestimate the 


half-widths of the R lines of C 2 H 2 broadened by N 2 . This 
demonstrates a known fact that the RB formalism predicts too 
large amounts of the half-width for two linear molecular sys- 
tems. Meanwhile, after considering the line coupling, the new 
formalism can significantly reduce calculated half-widths by 
amounts as large as 13%. Although, in comparison with mea- 
surements, there are still large gaps existing, but the refine- 
ment goes in the right direction. This conclusion has been 
found in our previous study on the isotropic Raman Q lines 
of N 2 broadened by N 2 and it is confirmed here again. 

Furthermore, in order to show the line coupling within 
the same branches plays a dominant role in determining its 
effects on calculated half-widths, we have also calculated the 
W matrix without including the line coupling between the P 
and R branches and present the calculated results in Fig. 7. As 
shown in the figure, the differences between these two results 
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R(j) 


FIG. 7. Comparison of calculated half-widths of the infrared R lines for C2H2 broadened by N2 derived from with and without taking into account effects from 
the line coupling. They are represented by □ (blue) and + (red), respectively. Meanwhile, some measured data 16 are plotted by x . In addition, theoretically 
calculated half-widths of the R lines by neglecting inter-branch coupling are plotted by (green). 


are small. More specifically, all the differences are less than 
2% and after j > 8, they become less 1%. Therefore, one can 
conclude that to ignore the line coupling between the P and R 
branches is a good approximation for this system. To isolate 
the R lines from the P lines would reduce the size of the S 2 
matrices by a half, but it can still yield good results. In the 
present study, because the size of S 2 (i.e., 42 x 42) remains 
well tractable, one does not need to rely on this simplification. 
But, it is worth to mention it here. 


IV. DISCUSSIONS AND CONCLUSIONS 

First of all, we would like to compare our new method 
with approaches used by others. It is well known that 
people have developed several formalisms to consider the 
line coupling for different molecular systems. Among them, 
Cherkasov has spent more than three decades to study this 
subject 20 and one of his papers appears recently. 24 In our 
previous study, 11 we have provided a detailed comparison 
of his method and ours. Meanwhile, there are two recent 
papers by Starikov. 26 ' 26 In the first paper, the author has 
introduced so called bi-resonance functions in the theory of 
collisional broadening of the spectral lines of molecules to 
consider the line interference. He has provided comprehen- 
sive lists of the forms with which these bi-resonance functions 
associated with the straight-line, the parabolic, and the “ex- 
act” trajectories can be modeled. In the second one, with the 
bi-resonance functions he has considered overlapping lines of 
the ammonia molecule broadened by argon and helium. 

It is worth mentioning that there is no essential difference 
between Starikov’s bi-resonance functions defined as summa- 
tions of products of resonance functions and their partners’ 
complex conjugates over the magnetic quantum numbers and 
those resonance products appearing in Cherkasov’s formulas. 


In addition, the two same arguments (i.e., the initial energy 
gap and the final energy gap of the two coupled lines) serve 
as the two modulation frequencies in both these resonance 
forms. Of course, Cherkasov only considers simpler potential 
and trajectory models and Starikov covers more models. With 
respect to a more profound problem resulting from how to di- 
agonalize matrices with a large size at each step in averaging 
over the states of the bath molecule, Starikov has completely 
followed the approach by Cherkasov. As explained in our pre- 
vious work, 11 by correctly applying the cumulant expansion, 
one is able to dramatically reduce the size of matrices to be di- 
agonalized. As a result, an extra approximation introduced by 
Cherkasov to separate the relaxation operators into one diag- 
onal part and another off-diagonal part becomes unnecessary. 

Furthermore, by fully exploiting the coordinate represen- 
tation, the 2D correlation functions can be introduced and by 
making variable changes these functions become even. Then, 
the latters’ 2D Fourier transforms are real and the complex- 
ity resulting from involving complex functions in calculating 
the off-diagonal matrix elements is completely obviated in 
advance. In addition, in comparison with others’ formalism 
given in terms of the bi-resonance functions and the products 
of resonance functions, contributions from one specified 2D 
Fourier transform ^l x l 2 (a>, co ') alone in our formalism corre- 
spond to a combination of contributions from a whole group 
of those resonance forms. The group categorized by the com- 
mon rotational symmetry labels of Li and L 2 consists of a lot 
of the resonance forms having different other indices, includ- 
ing those used to distinguish different R dependences. This 
implies that with our formalism, we are able to dramatically 
reduce the number of these functions. In fact, we believe that 
the number has been reduced to the minimum. Finally, the 
two modulation frequencies of these 2D Fourier transforms 
represent real intrinsic connections between two coupled 
lines. All these features uniquely exist in our new formalism. 
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Concerning with our numerically calculated half-widths 
of both the isotropic Raman Q lines and the infrared P and 
R lines for C 2 H 2 broadened by N 2 , the conclusion is that ef- 
fects from the line coupling are important. This conclusion 
is well expected because the C 2 H 2 molecule has a small ro- 
tational constant. In general, in comparison with that derived 
from the modified RB formalism, calculated half-widths with 
the line coupling are reduced by 12% — 14% and become closer 
to measurements. 

However, if one investigates the whole relaxation matrix 
W further, one needs to verify the detailed balance principle 
(which connects the elements of W' f ,l f and W'^' f ) and the 
sum rule 27 (which connects the diagonal and off-diagonal el- 
ements of W). For the isotropic Raman spectra of N 2 , these 
properties were recently studied and the limits of a semi- 
classical treatment were analyzed. 22 A similar study for the 
infrared spectra is under development. 

Finally, we would like to mention some procedures that 
have not been taken into account in the present study. Be- 
cause the current potential model does not contain any vi- 
brational dependence, we are not able to meaningfully cal- 
culate the imaginary term iSi which is a main component of 
the imaginary part of iS 1 + S 2 for vibration-rotational spec- 
tra. In addition, it is well known that the imaginary part of 
iS 1 + S 2 is much smaller than the real part. As a result, we 
have neglected this imaginary part in the present study. There 
are no any problems to include the imaginary part in our fur- 
ther calculations. More explicitly, we know how to derive the 
matrix elements of ImS2,outer,i and ImS2, 0 uter,f because we have 
an effective tool at handy and have had experiences to carry 
out these calculations. 29 In fact, by introducing the 2D causal 
correlation functions 30 and by carrying out the corresponding 
2D Fourier transforms, one is able to obtain the 2D Hilbert 
transforms with which one can easily calculate these imagi- 
nary matrix elements. Meanwhile, if the potential model con- 
tains vibrational dependences, the correlation functions in- 
troduced in our formalism would depend on the vibrational 
quantum numbers as well. This implies that one may need 
to have different sets of the correlation functions in calculat- 
ing the matrix elements for each of the S 2 , outer, i, S 2 , outer, f, and 
S 2 , middle terms. We expect that whether to distinguish these 
different sets is necessary or not would depend on cases of 
interest. 

As a next research step, we will consider systems involv- 
ing more complicated molecules. The target is to consider the 
line coupling for the H 2 O-N 2 system because a lot of values 
of the half-width and shift of H 2 O lines listed in HITRAN 
come from theoretically calculated results with the RB for- 
malism. For this important system, we expect that the line 
coupling would happen among lines belonging to different 
branches and effects on calculated half-widths and shifts from 
the line coupling could be significant. At present, after suc- 
cessfully considering the line coupling for all three branches 
of linear molecules, we stand on a more sound position to 
pursue this difficult problem. 
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